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0^ ' Three-nucleon forces have been considered to be one possibihty to 

0^ ' resolve the well known discrepancy between experimental values and the- 

oretical calculations of the nucleon analyzing power in low energy nucleon- 
deuteron scattering. In this paper, we investigate possible effects of two- 
pion exchange three-nucleon forces on the analyzing power and the differ- 
y , ential cross section. We found that the reason for different effects on the 

analyzing power by different three-nucleon forces found in previous calcu- 
lations is related to the existence of the contact term. Effects of some vari- 
ations of two-pion exchange three-nucleon forces are investigated. Also, 
an expression for the measure of the nucleon analyzing power with quartet 
P-wave phase shifts is presented. 



Abstract 



Differential cross sections for nucleon-deuteron elastic scattering have peaks 
at forward and backward scattering angles and a minimum at a cm. scattering 
angle of, e.g., 9 ~ 105° at E^^i^ — 3 MeV. Around the cross section minimum 
angle, some observables calculated with realistic nucleon-nucleon (NN) poten- 
tials are known to deviate systematically from experimental values [Q. The 
nucleon analyzing power Ay (9) for energies below « 30 MeV has exhibited a 
notable discrepancy |g, Q], which is referred to as the Ay (9) puzzle. E.g., in 
the neutron-deuteron (n-d) elastic scattering at E^^b = 3 MeV, experimental 
Ay (9) has a maximum value at ~ 105° Q, while theoretical calculations with 
modern realistic NN potentials |^, |[ 0] undershoot the value by about 30 %. 
The three-nucleon (3N) system has been considered as a good testing ground for 
the NN interaction. The discrepancy between the experimental and calculated 



Ay (9) may show that there is room for improvement of modern NN potentials. 
Actually, it was pointed out that changes in ■^Pj NN forces or the spin-orbit 
component of a potential cause a dramatic increase in Ay{9) & 0, 10, O. 



However constraint from NN observables made it difhcult to obtain reasonable 
changes in the NN potential to resolve the Ay{9) puzzle |l^, |l^, Q. 

Another possibility for resolving the Ay {9) puzzle is the introduction of a 
three-nucleon force (3NF) into the nuclear Hamiltonian. It is well known that 
most realistic NN forces underbind the triton, and a 3NF based on exchange 
of two pions among the three nucleons (27rE-3NF) can explain the needed at- 
traction. So far, several 27rE-3NF models have been proposed, among which 
the Tucson-Melbourne (TM) 3N F p|] and the Brazil (the earlier version, BR' 
[ [l6| , and the latter version, BR [|17[) 3NF have been used for 3N calculations. 
Although these 3NF models were made based on different ideas in constructing 
off-shell ttN scattering amplitudes which are important ingredients in 27rE-3NF, 
the resulting potentials have essentially the same form with slightly different 
parameters. It is reported that with introducing the TM-3NF or BR-3NF, the 
calculated Ay (9) decreases, which means that the discrepancy with the exper- 
imental value is enhanced ||l^, Q. On the other hand, the calculations with 
the BR'-3NF or another 3NF model, the Urbana (UR) 3NF, are reported to 
improve Ay (9) slightly [ po|, pl[ . The UR-3NF is based on the A-mediated two- 
pion exchange diagram [ p2[ , which is a part of diagrams included in TM-3NF 
and BR-3NF. The discrepancy of the effects on Ay (9) should arise from a struc- 
ture difference between TM/BR-3NF and BR'/UR-SNF. In this paper, we study 
effects of the 27rE-3NF on Ay (9) carefully and investigate the possibility for re- 
solving the Ay (9) puzzle with a 3NF. All calculations are performed at E^^^i, = 
3 MeV, where experimental data are available for the differential cross section 
[|3| and Ay {9) ig. The Argonne Vis model (AV18) |§ is used as the input NN 
potential throughout this paper. 

Our method for calculating the 3N continuum state is based on a natural 
extension of the bound state calculation |gj, |2^, |l9| , in which the Faddeev equa- 
tion is expressed as an integral equation in coordinate space. In the continuum 
state calculation, there appear additional singularities in the Faddeev integral 
kernel, which are absent in the bound state calculation: elastic singularity and 
three-body breakup singularity. The latter does not appear at energies below 
the three-body breakup threshold as in the present work. The former singular- 
ity can be easily treated by the usual subtraction method [E6|. In the present 
calculation, 3N partial wave states for which NN and 3N forces act, are re- 
stricted to those with total two-nucleon angular momenta j < 2. The total 3N 
angular momenta (J) is truncated at J = 19/2, while 3NF is switched off for 
3N states with J > 9/2. These truncating procedures are known to be valid for 
the low-energy {E^at — 3 MeV) n-d scattering. 

The two-pion exchange three nucleon potential has the following form in 



Table 1: Various parameters for the three- nucleon potentials, Eq. (|^), used in 
the present work. 



3NF a (^-1) b (/i-3) c ifi-^) d ifi-^) 



BR' 
BR 
BRa 



-1.05 
1.05 
0.00 



-2.29 
-2.29 
-1.49 



0.00 
1.05 
0.00 



-0.768 
-0.768 
-0.373 



momentum space: 

T^(q,q') 



(27r)6 



U 



2^ 2 12^ 2 ("'l ■ q)(cr2 • q ) 



(fi.f2){a + 6(q-q')+c(<z2+g'2)} 
+ (JT3 • f2 X Ti)(icr3 • q X q')d , 



(1) 



where q and q' are the momenta of the propagating pions, /i is the pion mass 
and F{q^) a form factor which is parameterized as the dipolc form with a cutoff 
mass A. The parameters, a, b, c, and d, for the BR'-3NF U& and the BR-3NF 
[ITI are shown in Table 0. Since the Brazil 3NF model is based on the effective 
Lagrangian approach, in which several diagrams are considered explicitly, we can 
separate out the 3NF component which results from A-mediated diagram. The 
parameters for this 3NF component, which should correspond to the UR-3NF, 
are shown in Table Q as BRa. 

The cutoff mass A is chosen so as to reproduce the triton binding energy. 
The value of 700 MeV is used for the BR'-3NF and the BR-3NF, and 800 MeV 
for the BRa-3NF. Hereafter these are designated as BRypp, BRjoo, and BRa^soo, 
respectively. 

In general, analyzing powers are defined as a difference between cross sections 
with different orientations of incoming particles normalized to unpolarized cross 
sections. Therefore, before discussing the n-d polarization observables, we make 
a comment on effects of a 3NF on the unpolarized n-d differential cross section 
(DCS). From calculations for various combinations of NN potentials and 3NF 
models, we found that the calculated values of the DCS around the minimum 
region {6 = 105°) have a correlation with those of the triton binding energy, ,63. 
This is shown in Fig. |l|, where we plot the calculated values of the DCS (105°) 
for E2g^i, = 3 MeV against the calculated B3. The n-d DCS consists of spin- 
doublet scattering, spin-quartet scattering, and their interference terms pj\ . 
The above correlation can be understood as a result of the well known relation 
between the doublet scattering length (^a) and B3: the Phillips plot. The 
S-wavc DCS at low-energy is proportional to l/(fc^ + 1/a^), where a is the 
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Figure 1: Calculated values of the n-d differential cross section at 6* = 105° 
for E2ci,i, = 3 MeV plotted against the calculated triton binding energy B3. 
Experimental value is taken from Ref . [E3| . 



scattering length and k is the momentum, which means that the DCS decreases 
if the scattering length a becomes smaller. In fact, the calculated value of the 
doublet scattering length, 1.35 fm for AV18 (-83 = 7.51 MeV), turns out to be 
0.68 fm for AVI8+BR700 (B3 = 8.36 MeV), while the quartet scattering length 
is unaffected by a 3NF. Thus the decrease of the DCS with the introduction of 
a 3NF should be attributed to reproducing the triton binding energy. In Fig. Q, 
we observe that with reproducing B3, the DCS (105°) gets closer to the central 
value of the experiment. However, due to a rather large error bar, it is not 
conclusive whether the decrease is favored. 

In Table |, calculated values of S3; ^^(105°) and DCS(105°) at E^ab = 3 
MeV are shown for AV18, AV18+BR'7oo, AVI8+BR700, and AV18+BRa,800, 
together with the corresponding experimental values. Ay (6) at 9 ^ 104.0° [Q 
and DCS at 6* = 103.9° ||. We observe slight decrease (increase) of ^^,(105°) 
for AVI8+BR700 (AV18+BRl,oo and AV18+BRa.80o) compared to AV18, which 
is consistent with previous calculations jig, p^ po| , [2l| . In Table g, results for 
a modified version of AV18 (Mod-AV18) |l^, in which factors 0.96, 0.98, and 
1.06 multiply the ^Pq, ^Pi, and ^^2 AV18 potentials, respectively, are also 
shown. The modification causes relatively large effects on NN analyzing power: 
overshooting of peak values of neutron-proton Ay (9) by about 30 % and 10 % 
for E2ni, — 3 MeV and 25 MeV, respectively, which has been strongly criticized 

The BR-3NF and the BRa-3NF give opposite Ay (6) effects. From Table g, 
we see that there is no term corresponding to the coefficients a and c in the 
BRa-3NF, which comes from the ttN S'-wave scattering amplitude. Thus it is 



Table 2: The results of the triton binding energy; the neutron analyzing power 
and the differential cross section at 9 = 105° for the n-d scattering at E^ab = 3 
MeV with AV18 + various 3NF and the modified AV18 (Mod-AV18). Experi- 
mental values are ^j;(104.0°) |] and ct(103.9°) ||. 





B3 (MeV) 


Ay{105°) (%) 


ct(105°) (mb/sr) 


Exp. 


8.48 


5.96 ± 0.13 


90.6 ± 2.7 


AVIS 


7.51 


4.29 


93.5 


AVlB+BR^oo 


8.44 


4.50 


89.2 


AVI8+BR700 


8.36 


3.62 


89.6 


AV18+BRa,800 


8.37 


4.43 


89.5 


Mod-AV18 


7.53 


5.11 


93.4 



Table 3: The results of the triton binding energy; the neutron analyzing power 
and the differential cross section a.t 9 — 105° for the n-d scattering at -E^^^ = 3 
MeV with AV18 + each term in BR-3NF. 





B: 


3 (MeV) 


Ay 


-(105°) (%) 


cr(105°) (mb/sr) 


AV18+BRa 




7.48 




4.17 


93.7 


AV18-hBR6 




8.25 




4.55 


90.0 


AV18+BRc 




7.56 




3.81 


93.3 


AV18+BRrf 




7.63 




4.14 


92.9 



interesting to see how each term in 27rE-3NF affects Ay (9). To see this, we 
calculate the n-d scattering at E^ab = 3 MeV taking into account each term 
corresponding to the parameter a, or b, or c, or d in BR700 in addition to AV18. 
Each potential is designated as BR^, BR;,, BRc, and BR^, respectively. The 
results are shown in Table ^. From Table ||, we see that Ay{105°) decreases 
for BRq, BRc, and BR^, but increases for BR^. Especially BR^ gives a large 
Ay{105°) effect. From this, it is concluded that the contribution from BR^ is 
larger than the one from the other terms in BR'-3NF, BRa-3NF, and UR-3NF to 
give a small increase in Ay (9), while the contribution from BRc is overwhelming 
in lowering Ay (9) in BR-3NF and TM-3NF. The BRc includes the so-called 
contact term, which was argued to be excluded to avoid an odd behavior of 
the 3NF at short range |0 , or from a viewpoint of chiral constraints pa] . It is 
remarked that BR'-3NF is obtained from BR-3NF with a prescription to remove 
the contact term: replacing the coefficient a by a — 2/x^c and setting c to be 
zero. Therefore we may express that the different Ay (9) effect of BR-3NF from 
that of BR'-3NF arises from the existence of the contact term. 



In the AV18+BR'7oo (AV18+BRa,8oo) calculation, the DCS(105°) decreases 
by 5 % (4 %) compared to the AV18 calculation, while Ay{105°) increases by 5 % 
(3 %). On the other hand, in the Mod- AVIS calculation, ^j^(105°) increases with 
little change in the DCS(105°). Thus there is an essential difference between 
effects on Ay (9) from the 27rE-3NF and from the modification of the "^Pj NN 
force. The former is an increase in Ay (9) simply due to the decrease of the DCS 
due to the effect of reproducing the triton binding energy. 

In Table ||, we observe that each term in BR-3NF gives quite different effects 
in Ay (6). Next, we investigate each effect of the four terms in the 27rE-3NF. 
To do so, we introduced only the a {b, c, (i)-term as a 3NF by varying the 
coefficient of a {b, c, d) to reproduce the triton binding energy. These 3NF 
models are designated as Wa, Wb, Wc, and Wd- The results are shown in Table 
Q together with the values of the coefficients. Although the change of sign in 
a and c, and of the magnitude in a compared to the original values in Table yj 
may be unnatural, these 3NF models might be useful as phenomenological ones 
which reproduce the triton binding energy within a restricted functional form. 
A variety of Ay (9) effects are observed from these 3NF models: a large increase 
due to Wa', a small increase due to Wt and Wc; a relatively large decrease due 
to Wd, besides the decrease of in the DCS(105°) due to the binding energy 
effect. It is remarkable that Wa seems to reproduce the experimental value of 
Ay(105°) quite well. However, it turns out that the deuteron tensor analyzing 
powers are modified improperly by Wa at the same time. In Fig. 0, we plot 
Ay (9) at E2^^ = 3 McV (a) and T2oi9) at E'f^^ = 6 MeV (b) calculated with 
AVIS (solid lines) and AV18+Wa (dashed lines). We see that the experimental 
data for Ay {9) are well reproduced with the introduction of Wa, and T2q{9) is 
significantly modified. Although there is no T2o{9) data for n-d scattering, recent 
precise measurements of tensor analyzing powers for proton-deuteron scattering 
are reported to be well reproduced by calculations without a 3NF [g9[. Thus 
such distortion of T2o{9) for n-d scattering may produce another "puzzle". 

It is found that the H/a-3NF gives a different effect on n-d polarization 
observables than other 3NF models. We remark that this difference can be seen 
in the n-d quartet P-wave phase-shifts: (54p , Stp^,^ , and Sip^ , to which Ay{9) 
is known to be sensitive. The relation between the n-d phase shifts and Ay (9) 
is quite complicated, but as derived in the Appendix, a combination 

- 4M4p^^^ - 5Af4p3^^ -I- 9M4p^^^ (2) 

appears in an expression for Ay (9), where Mip^ = exp {iSip^) sin {Si pj). For 
small phase shift differences, Eq. (0) is proportional to 

4A3/2_i/2 + 9A5/2_3/2, (3) 

where Aj_j/ = 5ipj — Sip ,. Eq. (0) is a convenient expression for Ay{9) being 
consistent with results of three- nucleon phase shift analysis pQ]. In Table 0, 
we list the calculated values of (54p A3/2-1/21 and A5/2-3/2 for some models 



Table 4: The results of the triton binding energy; the neutron analyzing power 
and the differential cross section ai 9 = 105° for the n-d scattering at E^ab = 3 
MeV with AV18 + Wa, Wb, W^ and Wd- 

Bs (MeV) Ay{l05°) (%) cr(105°) (mb/sr) 



AV18 + W^a 


8.49 


5.93 


89.0 


(a = -14.4/i-i) 








AV18+T4^b 


8.50 


4.64 


89.0 


{b = -2.90Ai"3) 








AY18+Wc 


8.49 


5.13 


88.8 


(c = -1.25^-3) 








AVlS+P^rf 


8.50 


3.65 


88.9 


(d = -3.10^-3) 









Table 5: Phase shift for the n-d *Pi/2 state and the differences A3/2-1/2 and 
A5/2-3/2J which are defined in the text, at E2ab ~ '^ MeV. 





^^PU2 


A3/2-1/2 


A5/2-3/2 


AV18 


24.2 


1.9 


0.1 


AVI8+BR700 


24.5 


1.9 


-0.2 


AV18+BR^oo 


24.6 


1.4 


0.4 


AV18+BRa,800 


24.5 


1.6 


0.3 


Mod-AV18 


24.0 


2.2 


0.2 


AV18+M^a 


24.9 


0.2 


1.4 



presented in this work. From Table ra we see that A5/2-3/2 ~ for most cases 
except for Wa, for which A3/2-1/2 ^ 0. Therefore, Ay{9) is proportional to 
9A5/2_3/2 (4A3/2_i/2) for Wa (the other 3NF models). The difference of the 
factors, 9 and 4, explains the reason why Wa gives a large increase in Ay (6) in 
spite of the same order of the phase shift differences, A5/2-3/2 and A3/2_i/2- 
However, the difference seems to affect incorrectly the deuteron tensor analyzing 
powers. 

In summary, we have studied the effects of the 27rE-3NF, and its variations, 
on some observables for n-d elastic scattering at low energy. We found that a 
contact term included in the 27rE-3NF gives a rather large Ay (9) effect. This 
is the reason why effects on Ay (9) by BR/TM-3NF are different from those 
by BR'/UR-3NF, which does not include the contact term. Ay (9) increases by 
about 5% with a 27rE-3NF model in which the contact term is eliminated. How- 
ever, this increase is essentially the result of a decrease in the differential cross 
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Figure 2: Ay{e) at £;£„(, = 3 MeV (a) and T2o(6') at Ef^^ = 6 MeV (b) calculated 
with AV18 (solid lines) and AV18+VFa (dashed lines). Experimental data of 
Ay (9) are taken from Ref. Q. 



section caused by reproducing the triton binding energy. This contrasts with 
the increase of Ay{6) by the modification of the '^Pj NN force, which is caused 
by a variation in spin-dependent cross sections. We found a phenomenological 
3NF model which reproduces both of the triton binding energy and Ay{6). This 
3NF originates from ttN S-wave scattering in the intermediate state with the 
strength parameter adjusted to reproduce the triton binding energy. But this 
3NF destroys the good fit of the tensor analyzing power at the same time. Since 
forces arising from the exchange of pious should have a tensor character, it seems 
natural that such forces affect not only the spin vector observables but also the 
spin tensor observables. A 3NF involving any mechanism other than 27rE, which 
might have a character of a spin-orbit forces as suggested from the modification 
of the ^Pj NN force, should be examined to resolve the Ay {6) puzzle. 



Appendix 

In the spherical base, Ay (9) is given by 

i{e)Ay{e) = ti{e) {T+^{e) + T^^io)) /V2 (4) 

with !{&) = Tr (MM'') , and I{0)T^{9) = Tr (Mt^M'') , where M is a transition 
matrix and r^ is a nucleon rank-1 spin operator, whose matrix elements in the 
channel-spin representation are 

<sv\T^\s'iy' > = (-1)2^+5-''' V2ss' 

X (ss'-z/j/'|1-k)| ^^2 1/2 1 }• (5) 



Here, n — \/2n + 1, and s is the channel-spin. 

With partial-wave amplitudes, M'^^^,^, , the transition matrix elements, Msus'u' , 
are given as Q 

X (s'£V'O|J^')M,i,,,rr(0,O). (6) 

Here, we apply some assumptions. 

A. We assume off-diagonal matrix elements of the partial wave amplitude to 
vanish: 

Mfis'i' = Ss^s'Se.e'M2s+ii,j (7) 

B. Since we are interested in '^Pj waves, we consider contribution only from 
s = 3/2. Then we have: 

.3 3 , „,^ , r 3/2 3/2 1 \ ,^, 



with 



A'U,,{e) = ^ l(^^z.m,|Ji.')(^^i^'0|Ji^') 



xM.e,,Yp{e,0) (9) 



J' .£' .m'^ 

xM4,r™'*(0,O) (10) 



C. It is remarked that the shape of I{9)Ay{9) for the n-d scattering is roughly 
given by sin0. This ^-dependence arises when {£,me,£\m[) = (1, ±1,0,0), or 
(0,0, 1,±1). For these cases, after evaluating the summation over z/, i'' and mg 
in Eq. (0), we obtain 

i{e)T^{e) ex KY{^{e,o)M:s^^^J2^'^'p.' 

The summation in Eq. (p^) is proportional to Eq. (0) 
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